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Abstract 

We show that confluent Darboux-Crum transformations with emergent Jordan states are 
an effective tool for the design of optical systems governed by the Helmholtz equation under 
the paraxial approximation. The construction of generic, asymptotically real and periodic, PT- 
symmetric systems with local complex periodicity defects is discussed in detail. We show how 
the decay rate of the defect is related with the energy of the bound state trapped by the defect. 
In particular, the bound states in the continuum are confined by the periodicity defects with 
power law decay. We show that these defects possess complete invisibility; the wave functions 
of the system coincide asymptotically with the wave functions of the undistorted setting. The 
general results are illustrated with explicit examples of reflectionless models and systems with 
one spectral gap. We show that the spectral properties of the studied models are reflected by 
Lax-Novikov-type integrals of motion and associated supersymmetric structures of bosonized 
and exotic nature. 


1 Introduction 

Propagation of light beams can mimic evolution of matter waves in quantum systems. This follows 
from the formal coincidence of the Helmholtz equation for monochromatic light, propagating in 
magnetization free medium, and the stationary Schrodinger equation in quantum mechanics [ 1 ]. In 
this manner, behavior of quantum systems can be simulated by optical systems, and, vice versa, 
concepts common in optics can find their way to quantum settings. 

Nowadays, there can be constructed optical materials where intensity of light is subject to a 
controlled gain and loss m , 0, 0. Unusual optical properties of these systems are reflected by 
refractive index that acquires complex values. The perfect match between gain and loss prevents 
the light from being exponentially dimmed or brightened and it is reflected by invariance of the 
refractive index with respect to combined transformations of space-inversion (P) and time-reversal 
(T). The light passing through such a material can exhibit some remarkable properties. Among 
them, let us mention power oscillations or non-reciprocity of beam propagation that are caused by 
non-orthogonality of Bloch states 0, m, m, 0 : m, violation of Friedel’s law of Bragg scattering 
[10] . unidirectional invisibility nn, m, im m or invisibility of defects in the periodic structure 
of optical crystals m, uni- 


i 


The peculiar features of the PT-symmetric optical systems can be captured by the Helmholtz 
(Schrodinger) equation with a non-hermitian Hamiltonian. Such Hamiltonian operators have been 
studied extensively in quantum mechanics for last two decades, see m m, ns] for review. It was 
understood that reality of their spectra occurs due to existence of an antilinear integral of motion 
that was identified with the PT operator in most cases. The lack of hermiticity of the energy 
operator challenges probability interpretation for a quantum system; the scalar product has to be 
redefined [20]-[22] and its explicit form is usually exceedingly difficult to find. 

As the solutions of the Helmholtz equation correspond to a purely classical object (amplitude of 
the electric field), the hassle of the probability interpretation is avoided. In fact, the optical systems 
with balanced gain and loss open the door for realization of interesting phenomena predicted by 
non-Hermitian (PT-symmetric) quantum mechanics. One of them corresponds to the transition 
between exact and spontaneously broken PT symmetry, where the behavior of a physical system 
at the vicinity of the exceptional spectral points was of the main interest m , EDi. m- 

The Helmholtz equation for linearly polarized light acquires the form 

+ d l + V(x,z) = 0, (1.1) 

where oj is the frequency of monochromatic beam, and T is electric field. The refractive index 
n(x) = no + 5n(x), where no > 0 and |<5n(x)| << no, is PT symmetric, i.e. 5n(x) = 5n(—x)*. 
In dependence of its sign, the imaginary part of 5n(x ) represents loss or gain. In the paraxial 

. tiquj 

approximation, we set ty(x,z) = e l ~~ z ij}{x, z), where i/j(x,z) is an envelope function slowly varying 
in 2 . Neglecting the term with d 2 'tjj, the Helmholtz equation reduces then to 

id z ip(x, z ) = [~d 2 x + V (s)]^(x, z), (1.2) 

where z denotes a depth of propagation of light in the crystal while x is a transverse coordinate. The 
potential V(x) is proportional to 8n(x), V(x) ~ 5n(x). The form of physically acceptable solutions 
of (II.211 depends on characteristics of the considered system; they can be square integrable in x, 
or quasi-periodic in x when V{x) is a periodic function. It is worth stressing that contrary to the 
PT-symmetric quantum mechanic^], the quantity V(z) = J'^° \ip(x, z)\ 2 dx is physically relevant 
as it denotes the power of the beam if>(x,z) when measured in the depth 2 within the crystal. 
Due to the complex refractive index, the power V(z) can be oscillating function in z as the beam 
propagates though the medium, see e.g. [6j. 

In the current article, we will have in mind primarily the framework of paraxial approximation 
described by (11.21) . suitable for the low frequencies where the term d 2 ip can be neglected. Let us 
notice that another situation is also of interest in the literature, where no paraxial approximation 
is employed eh, m- By separation of variables ^(x,z) = e lXz 'ip(x), the equation (11.111 reduces 
effectively to + n ^ uj 2 — A 2 ^ ip(x) = 0. Here, the inhomogeneities of refractive index are 
developed in the direction of propagation of light, contrary to (11.21) where the light propagates in 
normal direction to the inhomogeneities of n(x). Our results will be applicable in this scheme as 
well. 

Supersymmetric quantum mechanics (SUSYQM) provides powerful techniques for construction 
and analysis of quantum systems. It allows us to alter interaction terms of both exactly- and quasi- 
exactly solvable models without compromising their solvability properties. It relates the scattering 
characteristics such as reflection and transmission amplitudes of the two systems [26] . It can serve 

1 The Hamiltonian H with complex, PT symmetric potential ceases to be self-adjoint with respect to the usual 

scalar product (/, g) = f^° f* (x)g(x)dx. To recover physical relevance of the scalar product and the self-adjoint ness 
of the Hamiltonian, it is necessary to make a suitable redefinition of the former one, see e.g. [18], |25j . 


2 





to find integrals of motion of the new systems; in the case of finite-gap and reflectionless systems, the 
nontrivial (Lax-Novikov) integrals can be identified in straightforward manner [27l[28lf29 i 30]. The 
framework of SUSYQM can also be used very efficiently in the analysis of the soliton scattering, or 
to construct (multi-)soliton solutions of Korteweg-de Vries (KdV) and modified Korteweg-de Vries 
equations proceeding from the trivial solutions, see e.g. m , m- 

SUSYQM relies on the intertwining operators. These are frequently represented in terms of 
higher-order differential operators and identified as Darboux-Crum transformations known from 
the analysis of differential equations [33]. They are usually used to construct systems with altered 
number of discrete energie^E Requirement of regularity of a new system imposes some restrictions to 
the use and setup of Darboux-Crum transformations. Discrete energies that lie below the spectral 
threshold of the original system can be added by the first order Darboux transformations. To 
generate discrete energy levels within a bounded gap separating the spectral bands of an original, 
periodic system, Darboux-Crum transformations of higher order have to be employed, see [35]. 

Methods of quantum mechanics provide efficient tools for analysis of optical systems. They have 
been used in the optics for a few decades [36] - m- In recent years, popularity of the supersymmetric 
approach witnesses growing interest in the context of the PT-symmetric optical devices [43j-[48J. It 
allows to construct solvable models of optical crystals with reflectionless interfaces [43], the systems 
possessing unidirectional invisibility [36], or periodic crystals with invisible periodicity defects |45jFI . 

Developing further the pioneering ideas of von Neumann and Wigner [39], a supersymmetric 
generalization of the procedure for the construction of spherically symmetric scattering potentials 
that support bound states in the continuum was proposed in Refs. [50] . m General aspects of 
the technique, known in the literature as confluent Darboux(-Crum) transformation, have been 
analyzed in the series of papers [52] [53] [53] [55] [56], see also refs. [5X 58]. The extension of the 
generalized framework to the systems defined on the whole real axis leads to the construction of 
non-Hermitian systems. Those were considered in the context of optics e.g. in [59], [60] . 

Bound states in the continuum (BIC) can be observed experimentally in optical systems [6l] 
and their theoretical investigation in the context of PT -symmetric lattices was done, for instance, 
in refs. [62], [HE]- The bound states located on the threshold of the continuum spectrum of PT 
symmetric systems were considered in [35!. 

In the present work, by a systematic employment of SUSYQM, we will construct PT -symmetric 
optical systems where the refraction index is asymptotically real and periodic , however, there are 
localized complex periodicity defects. We will show analytically that the decay rate of the defects is 
related to the energy of the bound states induced by them; the defects fall off as x~ l or x~ 2 for the 
systems with BIC, whereas exponential decay takes place for bound states associated with discrete 
energies. 

The work is organized as follows. In the next section, we will present the main characteristics of 
the confluent Darboux-Crum transformations (also called double step Darboux transformations in 
the literature), the basic tool for construction of the systems with bound states in the continuum. 
We will explain its relation to the standard supersymmetric quantum mechanics and provide formu¬ 
las for the bound states and potential of the new system. We will review basic properties of periodic 
systems in section [3] where we focus on the PT-symmetry of Bloch-Floquet states. In section [3] we 
will focus on the PT-symmetric potentials that support bound states in the continuum. Results of 
the systematic analysis are illustrated on the examples in section [5] We discuss there, particularly, 
reflectionless systems that possess either visible or invisible defects. We consider one-gap systems 

2 They also can be used to eliminate some states in the continuous part of the spectrum, see 1341 . 
invisibility of defects in the crystal structure means that the asymptotic form of the wave packet outgoing from 
the defect coincides with the wave packet that would propagate through the undistorted crystal. There, the reflection 
coefficient vanishes and the change in the phase factor reduces to the trivial one (modulo period of the potential), 
see e.g. [45] . 
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that possess bound states in the continuum as well. In section [6j we discuss some peculiar proper¬ 
ties of reflectionless and finite-gap systems. In particular, we focus on the integrals of motion that 
reflect spectral properties of the systems manifested in two kinds of associated super algebras. One 
of them corresponds to a hidden bosonized supersymmetry. It is based here on Lax-Novikov inte¬ 
grals mixed up with confluent Darboux-Crum transformations. The other one is an exotic hidden 
nonlinear supersymmetry based on the extended matrix Hamiltonian and containing the extended 
number of integrals of motion in comparison with a usual N = 2 supersymmetric structure. The 
last section is devoted to discussion and outlook. 


2 Confluent Darboux-Crum transformation and Jordan states 


For any quantum system we obtain here a partner by employing a confluent Darboux-Crum trans¬ 
formation, and observe the emergence of the Jordan states in the construction. Such states will 
appear later in the structure of the Lax-Novikov integrals controlling the invisibility of the PT- 
symmetric defects. 

Consider a quantum system given by the Schrodinger Hamiltonian 


d? 

H = -d^ + v(x) ‘ 

Let ipo{x) be a (physical or non-physical) eigenstate of eigenvalue Eq, 

(H - E 0 )^ 0 = 0 . 

We use it to introduce the first order differential operators 

= — A* := -T-i*. = -A - * . 

dx ipo dx ipo ipo dx dx ipo 


By definition, 


1 


Ai/j o = 0, 

A linear independent from i/jq eigenstate of H of the same eigenvalue Eq we take in a form 


A* 1 — = 0. 
Wo 


i>o\x) = 


Mx) ( r 

\JxQ 


ds 

^o( s ) 


fl G C , 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 


where xq € R is some fixed point. Notice that the Wronskian W(f,g) := fg' — f'g of i/jq and i^q 
is unit, W{^q,^q) = 1. In what follows by ip^(x) we shall denote a function associated with a 
function ip(x) according to the rule (12.51) . 

The application of A and to V’o 1 ’ and (1 /-00) b 5 respectively, gives 


Ai/> o b = 1 


V’o' 


y 


= —00- 


( 2 . 6 ) 


By virtue of Eqs. (12.41) and (12.6|) . the first order operators (12.31) factorize the shifted Hamiltonian, 

H-E 0 = A i A. (2.7) 

The alternative product of A and A$ gives the shifted Darboux (SUSY) partner Hamiltonian H , 


d? 

H — Eq = AA^ = — — ^ + V(x) , with V(x) = V(x) — 2 (logifjo(x)))' 


( 2 . 8 ) 
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The nodes of ip o give rise to singularities in V(x) not present in V(x). When V(x) is required to 
inherit the regularity of V(x), ipo has to be fixed as a nodeless function. The intertwining relations 

AH = HA, A^H = HA # , (2.9) 

follow then from (12.71) and (12.81) . The operator H — Eq annihilates the states 1/ipo and (l/ipof. 

In dependence on the nature of Eq and ipa(x), the Hamiltonians H and H are isospectral or 
almost isospectral up to eigenvalue Eq, which can be present in one system but missing in the other 
one. Operators A and A$ realize a mapping between two-dimensional eigenspaces of the second 
order differential operators H and H, 


Aip E = i’E, A^e = {E - Eq)iPe , (2.10) 

for any E ^ Eq and corresponding eigenstates, H^e = EtpE , HipE = EipE- 

Notice that the zero mode ipQ of H — Eq is mapped by A into the zero mode 1/tpo of H — Eq 
and, analogously, the zero mode (I/i/jq) of H — Eq is mapped by AJ into the zero mode i/jq of 
H — Eq. However, the pre-image of the zero mode ipQ of H — Eq with respect to the action of A^ 
is not a zero mode of H — Eq. Similarly, the pre-image of the zero mode (l/ipo)' a of H — Eq under 
the action of A does not belong to the kernel of H — Eq. Instead, the indicated pre-images are the 
Jordan states of the corresponding operators H — Eq and H — Eq. Indeed, consider the equation 

A x = -(l/iPo) b . (2.11) 


Its solution is 
X(x) = -ip 0 (x) 


f XQ 


Vd (s' 


■(l/ipo(s)yds + a) =-jpo(x 


ipQ(r)dr 



ds + a 


( 2 . 12 ) 


where a S C and a € C are some constants. Acting on both sides of (12.111) by A^, one finds that 


(H — Eq)x = tpo , {H - Eofx = 0 . 


(2.13) 


This means that x( x ) is the Jordan state of H — Eq. Analogously, we find that the solution of 
equation A# A = ^, given by \{x) = -^y ^ f* g 7p 0 (s)ip 0 l, (s)ds + b^j , be C, satisfies {H - Eq) A = 

l/V’o > {H — Eq) 2 \ = 0. Without loss of generality, we set below a = 0 in (12.121) . 

Let us take E / Eq and consider the eigenstate equation 


(H-E)V(x\E) =0. 


(2.14) 


We suppose that E is sufficiently close to Eq and denote E — Eq = e. Assuming the analyticity 
of ^f(x\E) in E in vicinity of Eq, we look for solution of ()2.14j) close to ipo(x) in the form of the 
Taylor series in e, 


^(x; E) = ip 0 (x) + V — Xn(x ), 

n\ 

n =1 


. , d n V(x-,E). 
x n( x ) = - ^7 - \ E=E 0- 


Substituting this series into (12.141) . we find that Xn( x ) ar e defined recursively by 

{H - Eq)xi(x) = ipo(x), 

(H - Eo)xn{x) = nxn-\{x), 


(2.15) 
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n = 2,... . 


(2.16) 

(2.17) 














Therefore, Xi( x ) is just the Jordan state that appeared in Eq. (12. 13|) . Xi(x) = x( x )- The state 
Xn( x ) with n > 2 satisfies relations (H — Eo) n Xn( x ) = P^oix), (H — Eo) n+1 Xn( x ) = 0, and can be 
identified as a higher, n-th order Jordan state. 

Consider the system H e = --pjJs+V £ (x) generated from H by applying to the latter the Darboux- 
Crum transformation based on the eigenstates ipo(x) and 4 '(x]E), 


V e (x) = V(x)-2(logW(M x )^( x -,E)))" . 
Taking into account equation (12.1511 , in the limit e — > 0 we get 

V := lim V £ (x) = V(x) — 2 (log W(ipo(x), x(x)))" ■ 


£—>0 


Relation W(f,g) = / 2 ^(j) together with Eq. (|2.12l) gives us 

V(x) = V(x) - 2-^ log T{x) = V(x) - 4 J°f° + 2 




where 


.K + “ (« + <*) 

X(x) = [ ipQ(s)ds + a, 

Jo 


2 ’ 


(2.18) 

(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 


and from now on we set xq = 0. 

In correspondence with the Darboux-Crum construction, the Hamiltonian operators H , H and 




( 2 . 22 ) 


are almost isospectral. Since in the limit e —> 0, HT(x; Eq + e) = A(ipo + ex) = zAx = ^(1/V^o) , 
we have 

H-E 0 = B i B, H - E 0 = BB*, (2.23) 


where 


B = i] o-y-—, = — ~~j~do > with t/ 0 := (l/ipo? - 

ax rj o r/o ax 


(2.24) 


There holds Btjq = 0, -B?/ 0 b = 1 / t / o , B$(l/r)o) = 0, B^(l/r] 0 ) b = — 7 / 0 , and so, ker (H — Eq) = 
span {1/r/o, (1 /?7o) b } • From (12.231) it follows that 


BH = HB , B^H = HBK 

The H and H are intertwined then by the second order differential operators, 

(BA)H = H(BA ), (A^B^)H = H(A^B ^). 

The eigenstates ip and ip of H and H, Hip = Eip , Hip = Eip, for E ^ Eq satisfy 


ip = (BA)iP = (E 0 - E)ip + 
(A^B^)ip = (E - E 0 ) 2 ip. 


4’o 


Jo ^Pl(s)ds + a 


(ipfpQ - ipQlp'), 


(2.25) 

(2.26) 

(2.27) 

(2.28) 


Notice that in accordance with Eq. ()2.27f) and (12.241) . the eigenstate ipo of H of eigenvalue Eq 
is mapped into the corresponding eigenstate of H, 


ip 0 b 


__ ipo 

do JJipQ(s)ds + a 


(2.29) 
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To simplify notations, in what follows we re-denote —"00^ by d'o- 

Thus, the confluent double step Darboux transformation coincides with the second order 
Darboux-Crum transformation based on the eigenstate i/jo(x) of H, ( H — Eq)i/;q = 0 , and the 
associated with it Jordan state x( x ) that satisfies ( H — Eq)x{x) = i/jo(x), (H ~ Eo) 2 x{ x ) = 0 , see 

Eq. (EH). 

The described picture can be generalized further by applying the confluent Darboux-Crum 
transformations based on the eigenstate ipo and the associated Jordan states xi , • ■ ■ ,Xn of H, see 
(12.161) . ()2.17l) . A new Hamiltonian H n . H\ = H, will be intertwined then with H by the nth-order 
operators. Instead of analyzing such new systems following the line presented here, we just notice 
that a Wronskian formulation for confluent supersymmetric transformation chains was discussed 
e.g. in [ 55] . 


3 PT-symmetric systems 

We are going to consider PT-symmetric systems with either periodic or asymptotically periodic 
potential. Theoretical aspects of the complex and PT-symmetric periodic potentials have been a 
subject of intensive research, see e.g. [63, ESI ISL :SZi ESI [69]. We will focus on the construction 
of PT-symmetric systems by applying the confluent Darboux-Crum transformations to Hermitian 
periodic Schrodinger Hamiltonians. We suppose that the initial system is given by a real, regular, 
even T-periodic potential defined on the real line 


V(x + L) = V(x), V(-x) = V(x), V*(x) = V(x). 

The generic properties of the solutions of the equation 


“Jp + y(x) = 


E € 


(3.1) 


(3.2) 


are described then by Floquet’s theorem m- It tells that the two linearly independent solutions 
of (13.21) . which are neither periodic nor antiperiodic, can be written in terms of quasi-periodic 
functions Q. 

kL 

i/j±(x) = e ±lkx Uk(±x ), Uk(x + L) = Uk(x ), k € C, —^Z. (3-3) 

7r 

Here, we have taken into account that the spatial reflection operator P is the symmetry of H. We 
deal with the bounded (“stable”) solutions as long as k is purely real. Their energies belong to the 
interior of the allowed bands. The functions (13.31) with Im/c ^ 0 are unbounded (“unstable”) as 
they diverge exponentially when x goes towards plus or minus infinity. The corresponding energies 
are non-physical and form forbidden bands in the spectrum. When one of the solutions of (13.21) is 
periodic or antiperiodic, i.e. when ^ £Z, the second solution can either have the same periodicity 
or can fail to be quasi-periodic at allH. The former property corresponds to the solutions inside the 
allowed bands, while the latter characterizes the states at the edges of the allowed energy bands [fl 

4 We call a function f(x) quasi-periodic when it satisfies f(x + L) = cf(x) where c is a complex number. 

5 When (13.21) has a periodic or anti-periodic solution p(x), then the other solution y{x) satisfies y(x + 2 L) = 
y(x) + 7 p(x). The actual value of 7 depends on the concrete form of the potential in (13.21) (one can have 7 = 0), see 
EH for more details. 

6 The free particle system is a particular example where all physical solutions are periodic. Indeed, the Hamiltonian 

j2 q 

H = — is translation invariant and the solutions sin kx and cos kx are periodic with the period — for any real 
k ^ 0, while the constant solution with k = 0, at the edge of continuous spectrum has an arbitrary period. 
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In order to get PT-symmetric Hamiltonian H , we suppose that the eigenfunction ipo of the 
initial PT-symmetric Hamiltonian H is also an eigenstate of the operator PT, 

PTijj o = etp 0 , e 2 = 1, (H - E 0 )ip 0 = 0 , E 0 € R. (3.4) 

The requirement of PT-symmetry of H is equivalent to PT-^j log X{x)PT = -j^r\ogT{x) with 
Z(x) defined in (12.2111 . Taking into account (13.411 and (|2.21ll . this requirement can be met provided 
that a is purely imaginary, 

a = ia R , a R € R . (3-5) 

Besides, a R has to be (and can be) fixed in a way such that V is free of singularities, i.e. I(x) is 
nodeless. We suppose this to be the case from now on. 

Let us discuss in more detail the consequences of the requirement (13.411 for both stable and 
unstable solutions. 

First, let us suppose that ipo is a linear combination of the stable states g± with real quasi¬ 
momentum, 

g± = a± e ±lkx Uk(±x ), &;€R, k / n—, |a ± | = l. (3-6) 

In definition of g_, we used the fact that P is a symmetry of H. The phase factors a ^ are fixed 
such that there holds 

PTg± = g± > Pg±=9Tg±=g T . (3.7) 

This implies that any linear combination b\g + + b^g^ with real b\ , 62 is PT-symmetric. It is also 
nodeless provided that | 6 i| 7 ^ 162 h this follows from the fact that the real and imaginary parts 
of the latter linear combination are two independent solutions of (|3.2I) . and, hence, cannot vanish 
simultaneously. 

Secondly, we consider the case when ^ € Z. Let us take ipo as a linear combination of the 
states where at least one of them is periodic or antiperiodic and PT-symmetric. Let us denote it 
p(x). We have 

i/}q( x) = c + p(x) + C-q(x), p(x + 2T) = p(x ), (3.8) 

where the function q(x) can be written as q(x) = p(x) fg ^y^ds. As it was already noted, such 
states are associated either with the band-edge energies or they can correspond to specific en¬ 
ergy levels from the interior of the allowed energy bands as well. When PTp(x) = ep(x), then 
PTq(x ) = —eq(x). Hence, ipo is PT-symmetric provided that c+ € R and c_ € iM. up to a common 
multiplicative factor. 

Finally, let us consider i/jq as a linear combination of unstable states (13.3li with complex quasi¬ 
momentum, k = &2 ~ ik \, k\ > 0, k\, k 2 € R. One can show that ^2 = ^ where n is an integer. 
Indeed, let f(x) = e klx e lk2X Uk{x ) is an unbounded solution of (13.21) . Then f*(x) is also a solution. 
Considering asymptotic behavior of the two functions, we find that f(x) = cf*{x ) for a constant c. 
We get c = e ~ 2lk2X \ . Substituting x —>• x + L on the right-hand side, the equation is satisfied 

“fe \'*'J 

provided that &2 = The equation also implies that Uk(x) and u* k (x) differ just by a multiplicative 
constant, and, hence, Uk{x) can be fixed as a real function. We can write the fundamental solutions 
in terms of real functions 


/+ = e klX u k (x) , /- = P/+ = e klX u k (-x) , (3.9) 

that satisfy the following relations 


Tf± = f± , PTf± = . 


(3.10) 
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To get an eigenstate of PT, we have to take the linear combination 

F± = Cf+ ± C*/_ , PTF± = ±F± , C = ci + ic2, ci,C2 G K. (3-11) 

These states are nodeless if C 1 C 2 7 ^ 0. In such a case, the real and imaginary parts form the 
fundamental set of solutions and cannot vanish simultaneously at the same point. 


4 Properties of H 

In this section, we address the question of square integrability of To = —V’o 1 ’ and of the asymptotic 
behavior of V. We shall consider separately three situations, distinguished by the three qualita¬ 
tively different forms of V’o discussed in the previous section. In the first two cases, V’o will be 
associated with the energy Eq belonging to an allowed spectral band, whereas in the third case, Eq 
will be non-physical eigenvalue of FL belonging to the spectral gap. 


Let us start with the case where V’o is a linear combination of the quasi-periodic functions (13. 61) . 


ipo(x) = c + e lkx u k {x) + c_e lkx u k (-x ), 


(4.1) 


where c± are real constants. In order to show that To is square integrable, let us suppose that V’o 
is a periodic function with period £, i.e. ipo(x + £) = i/jq(x). This requirement is satisfied when the 
periods of e lkx and u k are commensurable. It is worth to emphasize that in some cases, £ may be 
reduced to the original period L of the Hamiltonian H. We have 

[ V’o (s)ds= fo e ^^(s)ds + fff] £ ^(s)ds=[f]Q 0 + fjl] £ ^(s)ds, (4.2) 

•J 0 


where Qo = f^'ipQ(s)ds and [|] is the integer part of f. We suppose here that the integral Qq is 
nonvanishing. Using (14.21) and boundedness of - 00 ? w e can see that To decays asymptotically as 

1/x, 


To = 


V’o 


V’o 


fo^o( s )ds + ia R [f ]Qo + J^Y^o( s ) ds + ia R 


for 


00 


(4.3) 


and, therefore, the eigenstate To is square integrable. 

For the second case, let us consider ij)Q as fixed in (13.81) , 


ipo = c + p(x) + c-q(x), ip 0 (x + 2nL) = ip 0 (x) + n^p(x), 7 / 0 , 


(4.4) 


where 7 and c± are constants. Notice that when 7 = 0, V’o reduces to the form of Eq. 63) and 
(14.21) and (14.31) apply in this case. When 7 / 0, we can make the same steps that led to (14.30 . 
yielding 



for |x| —>• 00 , (4.5) 


where Cj, j = 0,1,2,3, are constants. Hence, the function Tq is also square integrable in this case, 


To = 


V’o 


f*^(s)ds + ia R x 2 


00 , 


(4.6) 
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In contrast to (14. 3|) . it decays as 1/x 2 . 

In the last case, ipo is given as a linear combination of unbounded solutions ([3. lip . We fix 


ip 0 = F + = e klX u k (x) + e klX u k (-x ), 


(4.7) 


(the analysis for other choice of ipo would follow similar steps). To make the forthcoming compu¬ 
tation compact yet easy to follow, let us introduce the following temporal notations, 


x = 


x 

II 


, J{x) = e 2k ' x u k 


rx—xL 

t (x)u' k (x), G(x) = Uk{x)uk{-x), KL(x) = / J{s)ds , (4.8) 

Jo 


and then we can write 


rx rx rx 

/ (s)ds = / e 2kis u\(s)ds + 2 / u k {s)u k {—s)ds = 
Jo J -x Jo 





Ah 

1 

h 


^ ^2 e 2kinL ^ J(s)ds + e 2xLkl IC(x) - e~ 2xLkl JC{-x)j +2 £ Q(s)d, 
(q smh2k\xL + e 2xLkl IC(x) — e~ 2xLkl K,(—x)^j + 2 f G(s)ds, 


(4.9) 


where we integrated by parts in the second line and summed over n in the third line. Here, 
Q = (coth kL — 1) J ( s)ds is a number, whereas JC(x) = K{x + L) is a periodic function. We can 
see that (14.91) increases exponentially at large |x|. Thus, the eigenstate 'ho decays exponentially for 
x —>• ±oo, 


T 0 = 


2 kiu k (x) -kix 

u%(x)—2K(x)—Q 

2k 1 u k (-x) kix 

u k (-x)-2K.(-x)-Q 


X —>• OO , 
X —>• —OO , 


(4.10) 


and represents a quadratically integrable bound state of H. 

The Hamiltonian H constructed by fixing ipo as either (ED or ED has some remarkable 
properties. First, we can observe that the potential term V(x) coincides asymptotically with V(x). 
This can be seen easily from the relation 


H = H-4T 0 ^o + 2T^. (4.11) 

For (14.11) . where -00 is bounded and T 0 decays as 1/x, the function V — V disappears as 1/x. When 
we have (14.4p . ipo has linear-like behavior whereas To decays as 1/x 2 . Hence, the term V — V 
vanishes as 1/x 2 . The periodicity defects V — V are invisible ; apart of the bound state T 0 , an 
arbitrary eigenstate ip of H coincides asymptotically with an eigenstate ip of H where ip = BAip. 
The wave function does not acquire any phase shift when passing through the defect. 

These conclusions stem directly from the chain of equalities 

ip = BAip = (E 0 - E)ip + x — (ipip' 0 - ip 0 ip') (4.12) 

Jo Vo + « 

= (E 0 - E)ip+ ^o(ipipo - ipoip'), Hip = Eip , (4.13) 

where the second term in (I4.13P vanishes asymptotically (IV’ol an< 4 IV’I are bounded), implying 
ip\ x - H=oo = {Eq — E)ip. Let us notice that for ipo = const (in this case H corresponds to the free 
particle), the second term in (14.111) cancels out. 
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The invisibility of the periodicity defect in V does not take place when ipo is fixed as a linear 
combination of unbounded, exponentially growing states. The potential term V of H is asymptot¬ 
ically periodic, however, it does not coincide with V in general. We have 


V = 


V - 8kiUk(x) ( 

V — 8kiUk{—x) 


k 1 u k {x)+(u k {x))' _ (g) \ 

u k {x)-Q-2K.(x) (u^.(x)-Q-2K.(x))' 2 ) ’ 

( k 1 u k (-x)+(u k (-x)Y _ kiu'l(-x) _ 

V u k (—x)—Q—2IC(—x) (u^(-x)-Q-2K(-x)) 2 


x —>■ oo, 
x —>• —oo , 


(4.14) 


which is an L-periodic and real (uk is L-periodic and real, see () 3.10 [) and (14.7(0 function. Also, the 
function ip differs asymptotically from ip as the second term in ( 14.131) is not vanishing. Notice that 
when Uk is constant, then /C(x) = 0, Q = 0 and V coincides asymptotically with V. This is the 
case when H is the Hamiltonian of the free particle. 

For some specific choices of H and ipo, the intermediate, PT-symmetric Hamiltonian H can be 
identified with the original Hamiltonian H up to a constant displacement of the coordinate, 


H{x) = H - 2(V£M))' = H{x + c) , (4.15) 

where c preserves @ PT-symmetry of H(x). When this is the case, the term V’o/V’o has period L, 
that is only possible provided i/jq = e tkx Uk(x). The requirement that ip o is an eigenfunction of PT 
forces k to be real. Hence, we can get (14.151) for physical energy Eq and quasi-periodic ipo only. The 
Hamiltonian H is intertwined with H by B = see (12.241) . where r )o is a linear combination 

of e lkx Uk(x + c) and e~ lkx Uk(—x — c). To make it satisfy PTrj 0 = 770 , we take 

770 = e lkx {c\Uk{x + c) + c[PTuk(x + c)) + e~ lkx (c 2 Uk{—x — c) + c^PTuk(—x — c)). (4.16) 


The Hamiltonian H = H(x + c) — 2 (t7q/?7o) / is not necessarily L-periodic. However, it can have 
period l when the periodicities of e lkx and Uk(x) are commensurable. When this is not the case, V 
fails to be periodic at all. In the next section, we will discuss some examples where these conclusions 
will be illustrated explicitly. 


5 Examples 


5.1 Reflectionless Systems 


Let us analyze the systems generated by confluent Darboux-Crum (double-step Darboux-Jordan) 
transformation from the free particle, 



(5.1) 


We will discuss the cases where ipo is a linear combination of either stable or unstable solutions 
of the equation (H — Eq)iI> = 0. First, we fix ipQ as a linear combination of the unstable states 
corresponding to the energy Eq = — k^ < 0. Without loss of generality, we can choose V’o i n the 
following form, 

ip 0 = cosh(A:ox + ir), ko > 0, r G 1. (5-2) 


Then, according to Eqs. (12.201) . (12.221) and (|2.29l) . the explicit form of the Hamiltonian El and of 
its quadratically integrable bound state Tq = —ipo^ can be written as 




cosh 4 (&o£ + *t) ko sinh(2(fcox + ir)) 


SI 


Sh 


,T, cosh(fe 0 x + IT) ^ 

_ - Bh -’ 


'Hence, c is either real or its imaginary part coincides with imaginary period of H{x). 
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where 


X 1 

Sh = ioL + - + —— sinh(2(/c 0 x + ir)) 
2 4fco 


(5.4) 


and a = a R — sin 2r. The amplitude of both the potential V and the bound state 'I'o decays 
exponentially for large |x|. 

As a second case, we consider ip o to be a linear combination of stable solutions of (H — Eo)ip = 0 
with Eq = fcg > 0 by mixing the plane waves e lk ° x and e ~ ik ° x , 


ip o = cos(kox + ir), k 0 > 0, r € K . 


(5.5) 


The function ipo is even with respect to the action of PT. The explicit form of the Hamiltonian H 
and of its bound state 'Fq is El 




cos 4 (/co x + ir) ko sin(2(kox + ir)) 


Q 2 

tJrji 


+ 


S T 


= 


cos(kox + ir) 

Ifr 


X 1 

St = ia + - + -rr- sin(2(fc 0 x + ir)), 
2 4/cq 


(5.6) 

(5.7) 


where a = a R — sinh2r. In contrast with the previous case, the amplitude of both the potential 
and the bound state decays now as l/x for large |x|. 

Next, we choose ipo as a single stable quasi-periodic eigenfunction of H, 


ipo = e ikoX . 


(5.8) 


In this case, H can be written as 


d 2 2&q 

dx 2 sin 2 (kox + ir) ' 


T = 2 1 °S ( 1 + 2k o a R ), 


(5.9) 


which corresponds to the PT-regularized trigonometric Poschl-Teller potential [72jEI- For this choice 
of ipo, the intermediate Hamiltonian H reduces to H. Hence, H and H are effectively intertwined 
by the first-order Darboux transformation. It is worth noticing that the periodic Hamiltonian H 
has a non-degenerate energy level at E = Eq. Indeed, one of the solutions of (H — Eq)J = 0 given 
by (12.291) is a bounded periodic function ’Po = —^ ipo^, 


’Fo = 


1 

sin(kox + ir) 


(5.10) 


The second solution grows asymptotically as |x| and, hence, it fails to be stable quasi-periodic 
state. 

Finally, let us take ipo as a linear combination of the band-edge states corresponding to Eq = 0, 

i.e. 

ipo = bx + ic, (5-11) 

where b ^ 0 and c are real parameters. Notice that (15.111) can be obtained from (15.21) or (15.51) in 
the limit ko —>• 0. As we found in the preceding section, in this case the potential V has to decay 

8 Some analogous but Hermitian system with singularities on the real line was considered in m- 

9 For a wider class of such systems see also m- 
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asymptotically as l/x 2 . The actual form of V and of the bound (quadratically integrable) state To 
is given by 


~ a9 2 (bx + ic) ( , (bx + ic) 3 \ ~ bx + ic 

——\ 2i —sH’ *»=— 

Sc = icx (bx + ic) H---h ia R . 

o 


(5.12) 

(5.13) 


For c = 0, we get the potential whose real singular analog was discussed e.g. in m in the context 
of solutions of the KdV equation. For b = 0, we get a regularized two-body Calogero potential 


V = 


2 c 4 


(c 2 x — ia R ) 2 


(5.14) 


Let us mention that a similar system, the PT-regularized Calogero model with both the centrifugal 
barrier and confining x 2 potential term, was considered in m , m, where its spectral properties 
were discussed in detail. 

The fact that H is intertwined with the free particle by the confluent Darboux-Crum trans¬ 
formation implies that the the systems described by H are reflectionless (the known reflectionless 
systems are related with the free particle system by the Darboux-Crum transformation, see e.g. 
HZ]). As we discussed in the preceding section, the systems constructed by fixing -i/>o as m or 
(15.111) possess invisible potential barriers since there is no phase-shift of wave functions. 

When ipQ is a linear combination of unbounded states (|5.2I) . then the potential V is still re¬ 
flectionless. However, when we substitute (15.2j) and ([5.311 into (|4.12(1 . ()4.13[) . we get the following 
asymptotic form of the wave functions in continuous part of the spectrum, 


lim BAe ±ikx = (E 0 - E + 2 k 0 (k 0 =F ik)) e ±ikx = (fc 0 T ik) 2 e ±ikx , (5.15) 

£—>•00 

lim BAe ±ikx = (E 0 - E + 2 k 0 (k 0 ± ik)) e ±ikx = (fc 0 ± ik) 2 e ±ikx . (5.16) 

£—>■ — 00 

Hence, the functions acquire nontrivial phase shifts by the interaction with the potential, that 
makes the potential detectable in principle. 

To illustrate our results, let us consider propagation of a light beam in the optical crystal with 
the refractive index described by V. For the reflectionless systems discussed in this section, one 
can find explicit solutions of the Helmholtz equation (11.21) . They can be obtained with the use of 
the intertwining operator BA and the solution of the “time”-dependent free Schrodinger equation 
id z i/j(x, z) = —d 2 i/j(x, z) that has the following form 

( X-XQ-VQZ) 2 

7/;(-T, 7) = ——- _ p i^o{x-XQ- 1 voz) p 2(^+2 iz) _ (5.IT) 

ycr 2 + 2 iz 

For 2 = 0, it coincides with the Gaussian wave packet i/j(x,0) = e luJ o{^-xo) e -(x-x 0 ) 2 / 2 a 2 _ fj erej 
xq is the initial position of the wavepacket, a is the width parameter while ojq and vq = 2 wo are 
the analogues of the wave number and and group velocity, respectively. Then the solution of the 
equation (11.21) with V can be constructed as follows (notice that the intertwining operator BA is 
z—independent) 

id z 'ijj(x, z) = [~d 2 + V(x)]$(x, z), 4>(x, z) = BAtp(x, z). (5.18) 


The intensity of the wave packet \ijj(x,z)\ 2 is depicted in the Fig. [H where the absence of any 
reflections on the potential barrier is manifest. In order to see the contrast with a reflective potential 
barrier, we plot in Fig. [2] the intensity of the wave packet propagating in the potential given by 
only the real part of a complex reflectionless potential. 
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Figure 1: Reflectionless optical potentials. On the left: PT-symmetric potential V{x) = 
“ 2 i& lo S (fo ^o( s ) ds + ia R ) discussed in section o The solid blue (dashed red) line represents 
the real (imaginary) part. On the right: intensity \ijj(x,z)\ 2 of the light beam when propagating 
through the optical potential. There is a bound state of energy with (a) Eq = — k%, (b) Eq = 0, 
(c) Eq = kl The periodic system (d) has no bound states. The respective ipo and parameters for the 
construction of the potentials in each case are given as follows, (a) ij )o = cosh(fcox + ir), a R = 1.15, 
ko = 0.7, r = 1.4 (b) i/jq = bx + ic, a R = —2.8, b = —0.26, c = —0.5 (c) ipo = cos(kox + ir), a R = 4, 
ko = 0.7, r = 1 (d) V’o = e ik ° x , a R = 1, ko = 0.7. In ■*/>( x , z), we fix a = 2, xq = —50 and cuo = 7 t/2 
in all the cases, see (15.171) and (15.181) . 


z 



Figure 2: Reflective optical potential. Intensity \tjj(x,z )| 2 of the light beam when propagating 
through the real part of the potential, Re[R(x)], in the model (|5.12l) . The parameters are the same 
as in (b) of Fig. [TJ The two main reflections that appear about x = 0 correspond to those produced 
by each potential well, see (b), Fig. [T] 
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5.2 One-gap systems 

Let us take now the initial system described by the one-gap Lame Hamiltonian, 
d 2 

H = — —2 + V (x), V(x) = 2k 2 sn 2 (x, k) — k 2 = —2 dn 2 (x, n) + 1 + k' 2 


(5.19) 


The operator has period 2 K, where K = K(k) is the complete elliptic integral of the first kind, 
K(k) = Jq^ 2 ( 1 + k? sin 2 (/>) _ 1 // 2 d</>, k E (0,1) is the modular parameter, and k' = \/l — k 2 is the 
complementary modular parameter 0. The spectrum of the Schrodinger operator H consists of two 
bands. The finite (valence) band is formed by the energies E E [0, k /2 ], and the infinite (conduction) 
band stretches over the interval E E [l,oo). The functions 


T „ , x H (x ± a) , 

(*) = Q ex P ItxZ (a)] 

solve the eigenstate problem 

(x) = E(a )^± (x ), E(a) = dn 2 (a, k) . 


(5.20) 


(5.21) 


Here 0, H, and Z are Jacobi’s Theta, Eta and Zeta functions, while parameter a can take arbitrary 
complex values. The functions dn (x, k), as well as cn (x, k) and sn (x, k) that will be used later on, 
are Jacobi elliptic functions. 

The actual value of a determines whether the function is bounded or unbounded. Let us write 
a = /? + * 7 , /3 ,7 E R. The wave functions are unbounded (non-physical) when the corresponding 
energies belong to the lower infinite band, E E (—oo, 0), or to the finite spectral gap, E E ( k' 2 , 1). 
This happens provided /3 E (0 ,K), and 7 = K' or 7 = 0, respectively. The two bands of physical 
eigenvalues correspond to 7 = [0, K'], and f5 = K for a finite (valence), and j3 = 0 for infinite 
(conduction) bands. The properties of the functions (15.20|) with respect to the parity and complex 
conjugation were discussed in detail in [35] . and their behavior under the composite PT transfor¬ 
mation is given by 


PT\V± + P 

(x) 


(x) 

for /3 E (0, K), 

(5.22) 


(x) 

= 

for 

P E (0, K), 

(5.23) 

PT\ 

(x) 

= Tf + * 7 (x) 

for 

7 € [0, K'}, 

(5.24) 

pt^2 

(x) 

= —(x) 

for 

7 G [0,K'], 

(5.25) 


where ([5.221) and (15.23[) correspond to the lower semi-infinite forbidden band and the finite spectral 
gap, while Eqs. (|5.24[) and (15.251) correspond to the finite and infinite allowed bands, respectively. 

As discussed in the preceding section, we can construct operators H with bound states of 
arbitrary energy by taking tpo as an appropriate linear combination of (15.201) . In the current case, 
we deal with the situation where analytical calculation of V and i/jq gets exceedingly difficult in 
general due to the term fg i/iQ(s)ds whose analytical form is rather unreachable. We focus here 
to the special cases where if ;0 is 2iL-periodic (dn(x,/c)), or 4iL-periodic (cn(x,/c) and sn(x, k)), 
and corresponds to the band-edge energies. The band edge states and eigenvalues are given by the 
relations 


H dn (x, k) = 0, H cn (x, k) = k' 2 ■ cn (x, k) , H sn (x, k) = 1 ■ sn (x, k) . (5.26) 

10 Lame potential (15.1911 possesses also a second, purely imaginary period 2 iK', K’ = K(k'), that is behind its 
one-gap nature. 
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Taking ipo as one of the band-edge states, we get 


H = H + 


4k 2 sn (x, k) cn (x, k) dn (x,k) Fe 


Se 


+ 2 


c2 


where we abbreviated 


F e 


k 4 sn 4 (x, k) , 

( — x — in 2 a + E(x, k), 

^0 = 

sn (x, 

«) 

k 4 cn 4 (x, k) , 

Se = s —(1 — k 2 )x + in 2 a + E(x, k), 

^0 = 

cn (x, 

k) 

dn 4 (x, k) , 

y ia + E(x, k), 

^0 = 

dn (x, 

«) 


(5.27) 


(5.28) 


Here, E(x,/c) is the incomplete elliptic integral of the second kind, E(x,/c) = (( dn 2 (s,n)ds. The 
extra potential terms (in addition to initial V(x) from (15. 19f) ) vanish for large \x\ (due to asymptot¬ 
ically linear behavior of Se), so that V coincides asymptotically with V. In [35], it was found that 
real periodicity defects, associated with discrete energies in the forbidden bands, induce phase shifts 
of the wave functions. In the current case, a PT -symmetric defect associated with the band-edge 
energy does not alter the asymptotic form of the wave functions. 

Despite the explicit analytical form of Ft for generic choice of i/jq is unreachable, we can still make 
some interesting conclusions. Let us discuss briefly the properties of the intermediate Hamiltonian 
H when ipo = Tf^x). As it was discussed in [35], in this case Ft coincides with the original but 
displaced Hamiltonian Ft, 

H(x) = H{x + a + iK') — Eq , (5.29) 

where K' = K(k'). Notice that tt(x + a + iK') is PT symmetric whenever a = i'y or a = K + i'y 
for 7 = [0, K'], i.e. when T^(x) corresponds to the physical state. We can fix the solutions of 
(. tt — Eq)/ = 0 in the following form, 


if>± = ^(x + a + iK , )±PT^(x + a + iK'), PT^ ± = ±^±, (5.30) 


where e equals either + or —. By the construction, ( 15.301) are eigenfunctions of PT. They are also 
quasi-periodic and bounded. In order to construct the Hamiltonian Ft, we fix B in terms of rjo 
which is a linear combination of and i/j-. However, taking into account that 


PT^(x + a + iK') 


— (x T a + iK '), a = i'y , 
e ±M(i-.z(o)) f « ( x + a + iK >^ a = K + i'y , 


one can check that it is sufficient to consider only r/o = ^^.(x + a + iK'). This gives 


(5.31) 


H = Ft - 2(log m )" = H(x + a + iK') - 2(log Vo )" 


H(x + 2ry), 70 = ^“(a; + a + iK '), 

H(x), r]o = T“(x + a T iK'). 

(5.32) 


Clearly, tt has no bound state as I /770 is not quadratically integrable. Darboux transformation 
based on 70 produces from tt either the initial Hamiltonian H(x) or its complex shifted copy 
H(x + 2 * 7 ). 

Finally, let us make a brief comment on the unbounded states (|5.20[) that can be written as 
f± = e klx e lk2X Uk 0 (x), k\,k2 € M. In the one-gap system, the value of k2 distinguishes the valence 
and conduction bands; it is equal to 0 for the states from the lower forbidden band, whereas = 97 ^ 
for the finite gap, see [35] for the details. 
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6 Hidden and exotic supersymmetric structures associated with 
finite-gap system H , and Jordan states 

Spectrum of a generic periodic Hamiltonian possesses an infinite number of gaps and bands. How¬ 
ever, there is a class of systems whose spectra contain just a finite number of gaps. If we denote 
the number of gaps as n, the free-particle Hamiltonian (15.11) represents the system with n = 0, 
whereas the one-gap Lame Hamiltonian (15.191) qualifies for n = 1. The peculiarity of having a finite 
number of gaps in the spectrum is associated with the existence of an integral of motion known as 
Lax-Novikov operator £. This integral is represented by a (2 n + l)-order differential operator. 

In correspondence with the Burchnall-Chaundy theorem [751 1791 1801 [ 8 Tj . its square is a polyno¬ 
mial in H, 

2n+l 

C 2 = U(H-E m ), [H,C] = 0. (6.1) 

m =1 

Here, E m are the energies corresponding to the edges of the allowed energy bands. 

The Lax-Novikov integral reflects spectral degeneracy of the finite-gap system. It distinguishes 
two states corresponding to the doubly-degenerate energies in the interior of the bands. Particularly, 
in the one-gap case, the Bloch states (|5.20l) of the same energy are eigenstates of £ of eigenvalues 
differing in sign. For the free particle (15. 1 [) . the Lax-Novikov integral coincides with the momentum 
operator £ = — It annihilates a constant, which is the wave function corresponding to zero 
energy. In the case of the one-gap Lame system, the integral acquires a rather nontrivial form of 
the third order differential operator, that annihilates the band edge states sn (x,n), cn (x,k) and 
dn(x,/c), see [271 [28l [35] for its explicit form. 

The operator H constructed from a finite-gap Hamiltonian H possesses the integral of motion 
obtained by a dressing procedure, 

C = BA£A i B*, [H,£] = 0. (6.2) 

When £ is of order 2 n + 1, then £ has order 2 n + 5. Making use of the intertwining relations, we 
find that £ satisfies the relation 

2n+l 

£ 2 = (H- E 0 ) 4 U(H-E m ). (6.3) 

m= 1 

The operator £ inherits the properties of £. Bloch solutions of the equation (H — E)ip = 0 
corresponding to the double-degenerate energy level E are also eigenstates of £ of eigenvalues 
differing in sign. When E corresponds to one of the non-degenerate band-edge energies E m , m = 
1,..., 2n + 1, the corresponding physical states are annihilated by this operator. Let us clarify the 
point where £ and £ differ from each other. For the purpose we shall analyze the kernel of £. 

We suppose that E m ^ Eq for m = 1,..., 2re + 1. There are 2n+ 1 states ifj m in the kernel which 
are images of the states annihilated by £, 

i’m = BAlpm , £V’m = 0, 771 = 1, ..., 2n + 1 . (6.4) 

Let us define auxiliary functions £2 and £3 that satisfy ££2 = ipo and ££3 = x- Then we find the 
remaining four vectors from the kernel of £, 

/o = %‘, h = % l l (6-5) 

f2 = Vo 1 [ TTE [ V’o(r)£ 2 (r)dr£s, f 3 = r^ 1 [ / ^o(r)T, 3 (r)drds. (6.6) 

Jo V’o (s) Jo Jo Ms) Jo 
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(6.7) 


Using the explicit form (12.31) and (12.241) of A$ and B\ one can check that there holds 

B*h = o 1 , A*B*f 2 = S 2 , A*B*f z = S 3 ■ 

The functions f\, / 2 and / 3 are the Jordan states of H associated with Eq, 

(H - EoYfj = -, (H-E 0 y+ 1 f j = 0, j = 1,2,3. (6.8) 

% 

As it was shown in [27, 28], the finite-gap systems possess a hidden N = 2 superalgebra graded 
by the parity operator P. The two supercharges, anticommuting with P, are given by C and 
iPC. It was showed there that the superalgebraic structure is preserved by the Darboux-Crum 
transformation based on the band-edge states. In our current case, the intertwining operators BA 
and A^B^ are PT -symmetric, however, they do not need to have a definite parity with respect to 
the space inversion P. Hence, as neither P[ nor £ have definite parity, P is prevented from being 
a viable grading operator. Instead, we can consider the PT symmetry operator, which commutes 
both with H and C. Let us define C\ = iC and £ 2 = iPTC. These operators anticommute with 
PT, and generate the N = 2 supersymmetry of H, 

2n+l 

[C a ,H] = 0, {C a , C b } = -5 ab (H - E n ) , a, 6 = 1,2. (6.9) 

72—1 

Notice that the sign of the anticommutator is changed when compared to the Hermitian case 
discussed, for instance, in [27 , 28] , As the intertwining operators BA and A^B^ commute with PT, 
this superalgebraic structure can be identified easily in the system described by H as well, 

2?2+l 

[4,4 = 0, {4,4} = -Sab(H - Eo) 4 l[(H-E n ), a, b =1,2, (6.10) 

772—1 


where the supercharges are defined as C a = BAC a A^BK The grading operator of the superalgera 
(IfUOl) is PT. 

Let us make a brief comment on the finite-gap systems where the intermediate Hamiltonian 
satisfies H = H(x + c), see (15.81) or ([5.291) . In these cases, the Lax-Novikov operator associated 
with H reduces to (2 n + 3)-order differential operator C(x) = B{x)C{x + c)B^{x). It satisfies 

C 2 = Y\%+l{H-E m )(H-E 0 y. 

The scheme given by the relations (16.91) and (16.101) corresponds to a bosonized supersymmetry, 
where no fermionic degrees of freedom are present in the system [82l 183 . 8 i| . However, by using 
the confluent Darboux-Crum transformations, an extended, exotic supersymmetry can also be 
constructed. Let us define a matrix Hamiltonian, 


n 


H 0 \ 

0 H ) ’ 


( 6 . 11 ) 


which resembles the supersymmetric Hamiltonian of the standard supersymmetry, see [26]. Instead 
of having linear supercharges, we can construct a pair of supercharges of the second order based 
on the confluent Darboux-Crum transformations, by taking the Pauli matrix cr 3 , cr| = 1, as the 
grading operator, 


Qi 


0 A*B* \ 
BA 0 ) ’ 
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0.2 — icr 3 Qi . 


( 6 . 12 ) 










The supercharges commute with the Hamiltonian and anticommute with < 73 , resulting in the fol¬ 
lowing N = 2 non-linear superalgebra, 

[n,Qa]= 0 , {Qa,^} = 0, {Qa,Qfe} =2,5 ab (^- J Eo) 2 , 0,6 = 1 , 2 . (6.13) 

The hnite-gap nature of the Hamiltonians H and H guarantees the existence of integrals C and C. 
They can be used to define two bosonic supercharges for Ti, 


■&i = 


{H — Eq) 2 C 0 \ 

0 l ) ' 


= i(T 3-Sfi • 


(6.14) 


In this definition, the polynomial (H — Eq ) 4 was inserted to keep the same differential order 2n + 5 
of the diagonal elements (recall that £ is of order 2n + 1 whereas C is of order 2re + 5). The bosonic 
and non-linear property of the integrals «£? are summarized in the following superalgebra, 

2n+l 

[H,J? a \=0, [S? a ,* 3 \=0, {J? a ,J? b } = 26 ab (U-E 0 ) 4 '[[(n-Em), a, 6 = 1 , 2 . (6.15) 

m =1 

The full supersymmetric structure gets enlarged if we consider also a new set of fermionic integrals 
by multiplying the integrals (16.121) by the extended Lax-Novikov operator 2z?i, 


which satisfy 


Si — Qi&i, S 2 — iasQi^i , 

[n,S a ] = 0, {S a ,a 3 } = 0. 


(6.16) 


(6.17) 


Finally, the nonlinear superalgebra can be completed by means of the remaining (anti)commutation 
relations, 


2n+l 

iQ a ,J? 2 ] = -2 e ab S b , [S a ,S? 2 ] = - 2 e afe JJ (^ - E m ){U - E 0 ) 4 Q b , (6.18) 

m= 1 

{Q a ,S fe } = 2 S a Qb + 25 ab {E - Eo) 2 ^ . (6.19) 

Hence, the extended system (16.111) is characterized by exotic supersymmetric structure described 
by the four supercharges Q a and S a , and the two bosonic integrals 2 z? a , which together generate 
the nonlinear superalgebra (16.131) . (16.150 . f[6. 17j) . (|6.18D . (|6.19l) . 

7 Discussion 

We showed that the confluent Darboux-Crum (double-step Darboux-Jordan) transformation is 
an efficient tool for creating PT-symmetric systems that are asymptotically real and periodic, 
and have a periodicity defect disappearing for large |x|. We pointed out the difference between 
the confluent and usual second-order Darboux-Crum transformations. Whereas the intertwining 
operators of the standard one annihilate two (formal) eigenstates of the initial Hamiltonian, the 
intertwining operators of the generalized transformation annihilate, besides an eigenstate of H, also 
the associated Jordan state, see Eqs. (12.251) . (12.241) . (12.111) and (|2.13l) . 

The described confluent Darboux-Crum transformations were applied to a generic Hamiltonian 
with real even periodic potential (13.21) . We discussed the general aspects of the construction like 
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existence of (quadratically integrable) bound states and asymptotic behavior of the created systems. 
We showed that bound states in the continuous part of the spectrum are associated with invisibility 
of the periodicity defects. 

It was also showed that the decay rate of the periodicity defect is determined by the position 
of the bound-state energy in the spectrum; when it belongs to the interior of the energy band 
or it can be identified with the band-edge energy, amplitudes of the defects disappear as aT 1 or 
x~ 2 . When the bound state corresponds to discrete energy in the energy gap or lower forbidden 
band, exponential decay of the defect takes place. It would be interesting to verify whether this 
observation is of general validity, i.e. whether a periodicity defect of the 1/x or 1/x 2 decay induces 
a bound state in the energy continuum, whereas the defects with exponential decay are responsible 
for discrete energies. 

The application of the generic results was focused on the special class of periodic systems that 
possess finite number of gaps in their spectra. In the subsection 15.11 we considered reflectionless 
systems derived by the confluent Darboux-Crum transformation from the free-particle model. In 
particular, we constructed reflectionless PT-symmetric systems which are completely invisible. 

In the subsection 15.21 we discussed asymptotically periodic systems with periodicity defects 
derived from the one-gap system described by the Lame equation. It is worth noticing that in the 
recent work [35], similar systems with soliton defects described by a Hermitian Hamiltonian were 
constructed from the one-gap Lame system using standard Darboux-Crum transformations. There, 
an analysis was provided how to construct new Hamiltonians with periodicity defects, that induce 
bound states with discrete energy in the gap or in the lower forbidden band, below the valence band. 
The subsection 15.21 extends those results with the use of the generalized (confluent) Darboux-Crum 
transformations. It allowed us to construct systems with bound state energy of arbitrary value. In 
particular, we constructed explicit models where the bound state was associated with the band-edge 
energy. We also discussed the specific case where the intermediate Hamiltonian coincides with the 
original one up to a complex shift of coordinates. 

In section [6] we discussed a set of integrals of motion associated with finite-gap systems, which 
give rise to bosonized and exotic supersymmetries. We showed in detail how the superalgebraic 
structure of the integrals of motion differs from the standard case when the supercharges based on 
the confluent Darboux-Crum transformation are taken into account. 

Despite illustrating our results on the finite-gap systems, we would like to stress that the results 
of section [J] apply to a broad class of real even periodic potentials. 

The invisibility of periodicity defects discussed in this paper resembles the Klein tunneling. This 
phenomenon occurs in relativistic quantum mechanics where spin-1/2 particles can tunnel through 
strong electrostatic barriers. When the particles are massless, the barrier becomes invisible for the 
particles in the sense that it has no effect on their dynamics, independently on its actual form. 
Although Klein tunneling was not observed for elementary particles, it is manifested in carbon 
nanostructures where dynamics is governed by one- or two-dimensional Dirac equation. There, it 
can cause absence of backscattering of Dirac fermions on impurities in carbon nanotubes, see e.g. 
[85] and references therein. In our case, the invisibility is very model-sensitive. On the contrary to 
the Klein effect, even a slight change of the potential can make it visible as the beam or particles 
start to scatter off it. The physics behind these two phenomena is distinct; whereas Klein tunneling 
stems on existence of solutions corresponding to antiparticles, invisibility discussed in this paper is 
rather the result of a fine interference of transmitted and reflected waves on the particular potential 
barriers. 
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